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We apply the coherent potential approximation (CPA) to a simple model for disordered supercon- 
ductors with d-wave pairing. We demonstrate that whilst the effectiveness of an electronic Van Hove 
singularity to enhance the transition temperature T c is reduced by disorder it is not eliminated. In 
fact we give a qualitative account of changes in the T c vs. doping curve with increasing disorder 
and compare our results with experiments on the Yo.8Cao.2Ba2(Cui_ c Zn c )307_<s alloys. 
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According to a widely shared point of view the high transition temperature, T c , of the superconducting cuprates 
O | is the consequence of a Van Hove singularity, in their normal state electronic structure, enhancing an otherwise 
!_i ■ unexceptional, that is to say weak, effective electron-electron attraction pL Clearly, if correct, this scenario would 
Oh' imply that relatively modest tinkering with such conventional mechanisms of pairing as electron-phonon interaction or 
spin-fluctuation could suffice to solve the central problem of high temperature superconductivity. This, conservative, 
view is supported by two main observations. Firstly, most parameter free, first-principles, calculations find a Van 
Hove singularity near the Fermi energy cf in these materials Secondly, the experimentally observed rise and fall 

of T c with doping can be interpreted, very naturally, as due to ep passing a Van Hove singularity. Of course, there 
I | are many objections to this, so called, Van Hove Scenario However, in the light of the emerging consensus that 

the Cooper pairs in these superconductors have d-like internal symmetry M two of these objections are particularly 
direct and damaging. The first one is that in the case of d-wave pairing the Van Hove mechanism is not effective. As 
it happens, this has been adequately answered by Newns et al. |§j|. The other one concerns the particular sensitivity 
of (i-wave pairing to disorder and will be addressed in the present letter. 

In the case of s-wave superconductors the influence of scattering by defects, which do not break time reversal 
invariance, on T c is governed by the Anderson's Theorem |J and hence the disorder averaged density of states, 
n(e). Thus, the effect of disorder on the Van Hove Scenario is merely a question of smearing the structure in n(e) 
corresponding to the singularity. By contrast for d-wave superconductors there is no Anderson's Theorem and hence 
' »-~S the role played by disorder is dramatically different. Not surprisingly, even without the complications introduce d by 



Van Hove singularities the problem has been studied only recently and there are still many open questions [h 



In what follows, we generalize the Coherent Potential Approximation (CPA) for disordered superconductors Jl6| to 
the case of eZ-wave pairing. In general, the CPA is a significant improvement on beyond the self-consistent Born 
approximation (SB A) Jl7j], currently in use for disordered <i-wave superconductors. Elsewhere [|l8| we shall discuss 
in detail the formalism, differences between the CPA and the SBA, and various implications for experiments on high 
temperature superconductivity. However, in this letter we will focus our attention on investigating, for the first time, 
the viability of the Van Hove Scenario in d-wave superconductors in the presence of disorder. 

Our arguments will be based on extended, negative Uij Hubbard model on a lattice, whose sites are labelled by i and 
j, with random site energies e^. Moreover, we shall work in the Hartree-Fock-Gorkov approximation which implies 
that the usual one particle Green function matrix, in Nambu space, at the Matsubara 'frequency' e„ = ^(2rt + 1) 
^ ' satisfies: 
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(ie n - e l + n)5u + t a Aij \ f G n (l,j;e n ) Gi 2 (l, j; e n ) \ fl 0\ , , 

A£ {ie n + ei-fi)Su-ta J {G 21 (l,j;e n ) G 22 (l,j;e n ) J ~ ^ \ 1 ) ' {1) 



where the hopping integrals i,j and the pairing potentials Ay will be taken to be non zero only when the sites i 
and j are nearest neighbours, [i is the chemical potential, and we shall often refer to the Greens function matrix as 
G(i,j;e n ). The above equations are completed by the self-consistency condition: 

a ^ =^4^ e "" ,?Gi2(i ' j;e " ) ' (2) 

n 

where f/y is an attractive interaction energy between two electrons on nearest neighbour sites and 77 is a positive 
infinitesimal. To simplify matters we have assumed that the Hartree term can be absorbed into the hopping integral 
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Uj and dropped it from equation (1). As usual equations (1) and (2) are to be solved subject to the requirement on 
the chemical potential that 

n = lE c4£ "" G ii^^ e ») > ( 3 ) 

n 

where n is the number of electrons per unit cell. 

Whilst, later, we shall refer to experiments on the substitutional alloys Yo.sCao.2Ba2(Cui_ c Zn c )307_,5 |1£] we do 
not wish to be very specific about the physical nature of the point defects represented by the site energies £j. Indeed, 
we are content to provide a reliable analysis of the simplest possible non trivial model. Thus we take them to be 
independent random variables defined to have values and —^5 with equal probability, i, on every site. As might 
be expected we shall be interested in the average of G(i,j;e n ) over the above ensemble. To calculate G(i,j;e n ) we 
shall make use of the Coherent Potential Approximation which is well known to be the mean field theory of disorder 
for problems similar to the one at hand pQ |. 

To define a tractable problem of interest we shall assume that the sites form a square lattice. Then for Si = 
for all i, in the normal state, where Ay = 0, the spectrum features 4, well known, Van Hove singularities at e = 0, 
resulting in logarithmic divergence: n(e) « —ln(e). Below a certain temperature T c there is a solution with Ay 
and d-wave symmetry. As was pointed out by Newns et al. Q in this model when n changes the Van Hove Scenario 
obtains as readily as in models where the symmetry of superconducting state is s type. That is to say T c rises to 
a maximum at n = 1 and then falls as n increases from n — to n — 2 for a fixed interaction strength CZy. As 
mentioned above, the technical question we shall answer in this letter is what happens to this behaviour when the site 
energies are not zero but randomly ^8 or 

Clearly, the relevant, generic, consequence of disorder is the smearing of structure in the density of states n(e). 
To gauge the extent of this in our model we calculated n(e) using the standard CPA procedure [P0[ . The results, 
illustrating the smearing of the Van Hove singularities at e = 0, are shown in Fig. 1(a) for various values of the 
scattering strength S. In Fig. 1(b) we also show the corresponding self energy E(e) which in CPA depends only on 
e but not on the wave vector k. For our simple model S(e„), at the Matsubara frequency e„, satisfies the following 
CPA equation §X§ 

E(en) = (^-E(cn))5(i,»;e n )(^ + S(en)) ■ (4) 

It will be useful to note that in the weak scattering limit equation (4) reduces to the Self-consistent Born Approximation 
(SBA) 

S SM y^G(M;e„) (5) 

and hence on the real energy axis ImE(e) ~ n{e). Interestingly, in the non-self-consistent Born approximation ImS(e) 
is logarithmically divergent at the van Hove singularity. In the disordered case one would expect such a singularity 
to be smeared out into a peak. Indeed the 5 = 0.6 curve in Fig. 1(b) is fully consistent with this expectation. The 
other curves show the significant deviation, the split Van Hove singularity for instance, between CPA and SBA. 

Let us now turn to the case where both superconductivity and disorder are present p6| . Although the full CPA 
program can be implemented for the problem defined by equations (1-3) and the specification of the site energy 
ensemble, it is convenient to make the approximation, valid when the coherence length £o is much larger then the 
lattice spacing, that the pairing potential A. y - does not fluctuate very much and replace it in equation (1) by its average 
value Aij |2l| . For conventional s-wave pairing this leads to the Anderson's Theorem Q which means that the only 
effect of disorder is to replace the density of states in the gap equation by its ensemble average n(e) plf . Thus, in 
the Van Hove Scenario it is clear that the smearing of the Van Hove singularity in Fig. 1(a) implies a weakening of 
the Van Hove enhancement of T c ~ exp[— 1/ (Uniep))]. The T c at optimal doping will thus be significantly reduced 
by disorder, even for s-wave pairing. 

As mentioned earlier the theory of disorder in d-w&ve superconductors turns out to be very different [[To] fl5| |. 
Nevertheless, the above simplification remains both valid and useful. In short we shall deploy the CPA to find 
G(i,j;e n ) for an averaged pairing potential Ay (old) and recalculate A y (new) using equation (2) with Gi2(i, j; e n ) 

replaced by G\i (i, j; e n ) repeating the process until convergence. In this way self-consistency on the average is ensured. 

To derive the basic CPA equations for disordered d-wave superconductors let us define the coherent Greens function 
G c {i,j; e„) = G(i,j; e„) by the equation: 
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j V A ii (* e « - M - S 22 (e„))(5 i; - t, 



where we did not introduce any off diagonal self-energies such as Ei2(e„) and £21 (en) because for the single site 
perturbations of our model they would turn out to be zero. The next step is to consider the scattering of the 
quasi-particles propagating according to G c (i,j; e n ) by the defects described by the potentials: 

±_( ±U \ / E u (6„) 



In a straightforward application of the CPA principles S(e n ) and therefore G c (i,j; e n ) is determined by the condition 
that these defects do not scatter on the average. After some algebra this leads to the condition, similar to the one in 
equation (4), that: 

Sn(e n ) = (^-Sn(e n ))G^ 1 (t l *;e n )(i* + Eii(e n )) . (8) 

We solved this equation numerically by iteration using the fact that G1i(i,j; e n ) = G] ) 1 (j, j; e n , fx) where G° is the 
solution of equation (1) with = for every i, = A,j and the renormalized frequencies e n and chemical potential 
fx are given by 

en = en-ImS u (e„) , fx = ix + Re£ u (e„) . (9) 

Note that equation (9) together with equation (8) constitutes a closed loop of relations which determine e n or alterna- 
tively Sn(e n ). Clearly, to solve this we need the Greens function matrix G (i, i; e n ). This was obtained numerically 
by a very efficient recursion method [^2f . 

To place above formulae in the context of previous work we note that in the weak scattering limit they reduce to 
those investigated in references |]l0| , |l3| -|l5| where Sn(e n ) was obtained by the SBA in equation (5). By contrast the 
CPA resums diagrams to all orders including some which have been investigated by Nersesyan et al. However, 

unlike them we find, analytically, that, for small 6, near the Fermi energy e w 0, n(e) w ^ 0j as m the 

earlier work of Gorkov and Kalugin ]10| . 

In Fig. 2 we display our results for the evolution of the density of states n(e) and the self-energy S(e) with disorder, 
as described by S, at T = 0. As expected the overall effect of disorder is to fill in the v-shaped dip in the density 
of states near ep. However, this process is rather intricate due to the competition between contributions in fc-space 
from near the point where A(fe) = and the Van Hove singularities at the saddle points. As a result, for small 6, 
the strong scattering at the Van Hove singularity, evident in Fig. 1(b), is suppressed in the superconducting state by 
the dip, but eventually, for large S, disorder wins. Namely £n(0) recovers and the gap is suppressed. Interestingly, 
the size of the 'gap' as measured by the distance between the two peaks in Fig. 2(a) remains roughly constant as the 
low energy dip is filled in. We have also performed finite temperature calculations relevant to tunnelling experiments 
|p3| , but we shall present these in a separate publication |l8| . 

Finally we investigate the dependence of T c on the band filling n. Our results for T c vs. n are displayed in 
Fig. 3(a). For orientation we also calculated T c as a function of n for an on site only, negative U Hubbard model 
with Uij = —\U\Sij and the same \U\ and hopping integrals t\j on the same lattice. The corresponding results are 
displayed in Fig. 3(b). Surprisingly, in spite of the delicate interaction between the Van Hove singularities and the 
very different behaviour of the quasiparticle spectra at low energies the two curves are largely similar. Hence, we 
may conclude that the Van Hove Scenario survives disorder almost as readily in d-wave superconductors as in the 
more conventional s-wave case. In support of the suggestion that this fact might have something to do with the 
superconducting cuprates we reproduce, in Fig. 4 the measurements of T c , as a function of oxygen deficiency for 
various concentrations of Zn, made on Yo.8Cao.2Ba2(Cui_ y Znj,)307_5 samples by Bernhard et al. |L9j. Clearly if the 
Zn concentration is regarded of a measure of disorder and oxygen deficiency as that of doping these curves are very 
similar to our theoretical results in Fig 3(a). Thus our calculations support the interpretation of the experimental 
data by Bernhard et al. |l9| as evidence for rf-wave pairing. 
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FIG. 1. Density of states N(E) (a) and self energies £(-E) (b) for a normal state with various disorder strengths 5. 



FIG. 2. Density of states N(E) (a) and self energies £(-E) (b) for a superconducting state with various disorder strengths S, 
calculated for \U\ = 3.5i and n = 1. 



FIG. 3. Critical temperature T c vs. band filling n for d (a) and s-wave (b) superconductors and a number of disorder 
strengths 8 (\U\ = 3.5t). 



FIG. 4. T c as a function of oxygen deficiency 5 for Yo.sCao.2Ba2(Cui_ !/ Zn !/ )307-i with y — (full squares), y — 0.02 
(triangles ), y = 0.04 (circles) and y — 0.06 (stars) [19]. 
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